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Abstract — This paper applies recent results on the robust 
stability of nonlinear quantum systems to the case of a Joseph- 
son junction in a resonant cavity. The Josephson junction is 
characterized by a Hamiltonian operator which contains a non- 
quadratic term involving a cosine function. This leads to a sector 
bounded nonlinearity which enables the previously developed 
theory to be applied to this system in order to analyze its 
stability. 

I. Introduction 

In recent years, a number of papers have considered the 
feedback control of systems whose dynamics are governed 
by the laws of quantum mechanics rather than classical 
mechanics; e.g., see [1] — [13]. In particular, the papers [10], 
[14] consider a framework of quantum systems defined in 
terms of a triple (S,L,H) where S is a scattering matrix, 
L is a vector of coupling operators and H is a Hamiltonian 
operator. 

The paper [15] considers the problem of absolute stability 
for a quantum system defined in terms of a triple (5, L, H) 
in which the quantum system Hamiltonian is decomposed as 
H = Hi + H-2 where Hi is a known nominal Hamiltonian 
and H-2 is a perturbation Hamiltonian, which is contained in a 
specified set of Hamiltonians W. In particular the paper [15] 
considers the case in which the nominal Hamiltonian Hi is 
a quadratic function of annihilation and creation operators 
and the coupling operator vector is a linear function of 
annihilation and creation operators. This case corresponds 
to a nominal linear quantum system; e.g., see [4], [5], [7], 
[8], [13]. Also, it is assumed that H 2 is contained in a set 
of non-quadratic perturbation Hamiltonians corresponding 
to a sector bound on the nonlinearity In this special case, 
[15] obtains a robust stability result in terms of a frequency 
domain condition. 

In this paper, we apply the result of [15] to a quantum 
system which consists of a Josephson junction in a resonant 
cavity as described in the paper [16]. This enables us to 
analyze the stability of this quantum system. In particular, 
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this enables us to choose suitable values for the coupling 
parameters in the system. For the parameter values chosen, 
we show that the quantum system is robustly mean square 
stable according to the definition of [15]. 

II. Quantum Systems 

The main result of the paper [15] considers open quantum 
systems defined by parameters (S, L, H) where H = Hi + 
H 2 ; e.g., see [10], [14]. The corresponding generator for this 
quantum system is given by 

g(x) = -i[x,H] + c(x) (i) 

where C{X) = \L^[X,L] + \[L\X)L. Here, [X,H] = 
XH — HX denotes the commutator between two operators 
and the notation t denotes the adjoint transpose of a vector 
of operators. Also, Hi is a self-adjoint operator on the under- 
lying Hilbert space referred to as the nominal Hamiltonian 
and H 2 is a self-adjoint operator on the underlying Hilbert 
space referred to as the perturbation Hamiltonian. The triple 
(S,L,H), along with the corresponding generators define 
the Heisenberg evolution X(t) of an operator X according 
to a quantum stochastic differential equation; e.g., see [14]. 

We now define a set of non-quadratic perturbation Hamil- 
tonians denoted W. The set W is defined in terms of the 
following power series (which is assumed to converge in the 
sense of the induced operator norm on the underlying Hilbert 
space) 
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Here Su = S$ k , H ke = C fe (C*) £ > an d C is a scalar operator 
on the underlying Hilbert space. Also, * denotes the adjoint 
of a scalar operator. It follows from this definition that 
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and thus H 2 is a self-adjoint operator. Note that it follows 
from the use of Wick ordering that the form (0 is the most 
general form for a perturbation Hamiltonian defined in terms 
of a single scalar operator £. 
Also, we let 
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and consider the sector bound condition 

f(c,crf(c,c)<^cc+Si 

T 

and the condition 

/"(CCT/"(C,C*)<<*2. 

Then we define the set W as follows: 

H 2 of the form (fJJ such that 
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W 



conditions (0 and (O are satisfied 



(7) 



Reference [15] also considers the case in which the 
nominal quantum system corresponds to a linear quantum 
system; e.g., see [4], [5], [7], [8], [13]. In this case, Hi is 
of the form 
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where M € C 
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is a Hermitian matrix of the form 



M 



Mi M 2 
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and Ml = Nl\, M 2 = M 2 . Here, the notation * denotes the 
vector of adjoint operators for a vector of operators. Also, 
# denotes denotes the complex conjugate of a matrix for a 
complex matrix. In addition, we assume L is of the form 



JVi N 2 



where JVj G C roX " and iV 2 G C* 
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Also, we write 

N 2 



(10) 
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As in [15], we consider a notion of robust mean square 
stability. 

Definition 1: An uncertain open quantum system defined 
by (S, L, H) where H = Hi + H 2 with Hi of the form 
©, H 2 G W, and L of the form (IToT > is said to be robustly 
mean square stable if for any H 2 G W, there exist constants 
ci > 0, c 2 > and C3 > such that 
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Here 



a(t) 
a*(t) 



vector of operators 
We define 

C = Eia + E 2 a* 

= \ Ei E 2 \ 



denotes the Heisenberg evolution of the 
; e.g., see [14]. 
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where Q is a scalar operator. Then, the following following 
strict bounded real condition provides a sufficient condition 



for the robust mean square stability of the nonlinear quantum 
system under consideration when H 2 G W: 
1) The matrix 

F = -iJM - ^JNUN is Hurwitz; (13) 
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where 
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This leads to the following theorem which is presented in 
[15]. 

Theorem 1: Consider an uncertain open quantum system 
defined by (S,L,H) such that H = Hi + H 2 where Hi 
is of the form dg), L is of the form ([T0j> and H 2 G W. 
Furthermore, assume that the strict bounded real condition 
([TJl i. ( TT4T > is satisfied. Then the uncertain quantum system is 
robustly mean square stable. 

In the next section, we will apply this theorem to analyze 
the stability of a nonlinear quantum system corresponding to 
a Josephson junction in a resonant cavity. 

III. The Josephson Junction in a Resonant Cavity 
System 

We consider a quantum system consisting of a small 
Josephson junction coupled to an electromagnetic resonant 
cavity. This system has been considered in the paper [16] 
where a Hamiltonian for the system is derived. The paper 
[16] considers the system as a closed quantum system defined 
purely in terms of a system Hamiltonian. We modify this de- 
scription of the system to consider an open quantum system 
model which interacts with external fields by introducing 
coupling operators for the system in order to apply the 
results of [15] given above. The coupling operator which 
is introduced is taken as a standard coupling operator for 
a resonant cavity coupled to a single field as well as a 
corresponding coupling operator for the Josephson junction 
coupled to an external heat bath. 

A Josephson junction consists of a thin insulating material 
between two superconducting layers. As in [16], we consider 
a Josephson junction in a resonant cavity. This is illustrated 
in Figure [TJ The following Hamiltonian for the Josephson 
junction system is derived in [16] 

U = i-[/'(n'-n') 2 -^cos^ + ^-(p' 2 + W V 2 ) 
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where n' =^,U' = U + hug 2 
ih. Here 



(15) 

i, and [q' 7 p'] = 
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= p + gV utin 1 q' = q 



where q and p are the position and momentum operators for 
the resonant cavity. Also, n is an operator which represents 
the difference between the number of Cooper pairs on the 
two superconducting islands which make up the junction. 
Furthermore, <f> is an operator which represents the phase 
difference across the junction. Note that compared to the 
expression for the Hamiltonian given in [16], we have 
normalized the expression (JT3J by dividing through by a 
factor of h in order to be consistent with the convention 
used in [15]. 

The quantities U, J', n, g, u> are physical constants 
associated with the junction and the resonant cavity. In 
particular, U is the charging energy of the Josephson junc- 
tion, J' is the Josephson energy of the junction, n is an 
experimental parameter relating to the gate voltage applied 
to the superconducting islands, g is a parameter related to the 
dimensions of the junction, and cj is an angular frequency 
related to the detuning of the cavity; see [16]. 

By a process of completion of squares, defining new 
variables n" — n' — n, p" = p' — g\fhH) and neglecting the 
constant terms, we can re-write the Hamiltonian as follows: 



\W P" *" 






" 




" q' 









p" 


u' 
h 







n" 










<t>' 



Now we define new operators a\ = {uq' + ip' r )/^/2huj, 
and a 2 = (</>' + in")/y/2 which satisfy the canonical 
commutation relations [ai,a*] = 1 and [0,2,012] = 1. Then, 
the Hamiltonian can be re-written in the form 



H= l - 
2 



where a 
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and M is a Hermitian matrix of the form 

L ° 2 J 

®. Note that in order to write the Hamiltonian in this form 
with the matrix M satisfying (O, it is necessary to apply 
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Fig. 1. Schematic diagram of a Josephson junction in a resonant cavity. 



these canonical commutation relations and neglect further 
constant terms. Also, the fact that the operators a\ and a 2 
commute is used to re-distribute terms within the matrix M. 
(Formally cii and a 2 are defined on different Hilbert spaces 
but following the standard convention in quantum mechanics 
we can extend each operator to the tensor product of the two 
Hilbert spaces and then these two extended operators will 
commute with each other.) 

We now modify the model of [16] by assuming that the 
cavity and Josephson modes are coupled to fields correspond- 
ing to coupling operators of the form 



L = 



This modification of the quantum system is necessary in 
order to obtain a damped quantum system whose stability 
can be established using the approach of [15] and is quite 
reasonable for an experimental system which will experience 
damping in both the electromagnetic cavity and in the 
Josephson junction circuit. 

In order to apply Theorem [TJ to analyze the stability of 
this quantum system, we rewrite ( fTSI l as 



H = H 1 +H 2 



where Hi = ^ [ a 

J' cL-2+a- 



M 



a* 



and H 2 



— j-cos( 2 J^ 2 )- That is, we have a quadratic nominal 
Hamiltonian and a non-quadratic perturbation Hamiltonian. 
Then, we define ( = a 2 /\/2 and 



— cos(c+n 

-cos( C 



C). 
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From this it follows that 

/'(CO*/'(CO<4^CC* 

and 7 = ^p. 

The numerical values of the constants 10, g, U, and J' 
are chosen as in [16] as ^ = 100GHz, g = 0.15, U = 
2.2087 x 10" 22 , J' = 3.6652 x 10 11 . Also, we calculate 
7/2 = 6.8209 x 10~ 13 . The parameters ki and k 2 will 
be chosen using Theorem [TJ in order guarantee the robust 
mean square stability of the quantum system. Indeed, for 
various values of ki and n 2 , we form the transfer function 
G KuK2 (s) = E#Y,{sI-Fy 1 D and calculate its 
norm. 

The Hoo norm of G Kl , K2 (s) was found to be virtually 
independent of n\ and so a physically reasonable value of 
n\ = 10 11 was chosen. With this value of Ki, a plot of 



h -2 \ 



versus k 2 is shown in Figure [2] 



From this plot we can see that stability can be guaranteed 
for k 2 > 2.2 x 10 12 . Hence, choosing a value of k 2 = 
2.5 x 10 12 , it follows that stability of Josephson junction 
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Fig. 2. Plot of ||G K1 , K2 (s)||oo versus K2- 

system can be guaranteed using TheoremQ] Indeed, with this 
value of k-2, we calculate the matrix F = —iJM—^JN^JN 
and find its eigenvalues to be -5.0000 x 10 10 ±3.3507 x I0 3 i 
and -1.2500 x 10 12 ± 1.4842 x W 3 i which implies that 
the matrix F is Hurwitz. Also, a magnitude Bode plot 
of the corresponding transfer function G KlyK2 (s) is shown 
in Figure [3] below which implies that ||G Kl!K2 (s)|| OC; = 
5.5554 x 10~ 13 < 7/2 = 6.8209 x 10" 13 . Hence, using 
Theorem[T] we conclude that the quantum system is robustly 
mean square stable. 
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Fig. 3. Magnitude Bode plot of G KljK2 (s). 

IV. Conclusions 

We have applied a previous result on the robust stability 
of nonlinear quantum systems to a quantum system arising 
from a Josephson junction coupled to an electromagnetic 
resonant cavity. This system involved a cosine function in 
the system Hamiltonian which was a suitable non-quadratic 
Hamiltonian leading to a quantum system with a sector 
bounded nonlinearity. For specific numerical values of the 



system parameters, it was shown that the robust stability 
result can be used to choose coupling parameters for the 
Josephson junction system in order to guarantee robust mean 
square stability. 
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